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ABSTRACT 

Orbits of known extrasolar planets that are located outside the tidal circular- 
ization regions of their parent stars are often substantially eccentric. By contrast, 
planetary orbits in our Solar System are approximately circular, reflecting planet 
formation within a nearly axisymmetric, circumsolar disk. We propose that or- 
bital eccentricities may be generated by divergent orbital migration of two planets 
in a viscously accreting circumstellar disk. The migration is divergent in the sense 
that the ratio of the orbital period of the outer planet to that of the inner planet 
grows. As the period ratio diverges, the planets traverse, but are not captured 
into, a series of mean-motion resonances that amplify their orbital eccentricities 
in rough inverse proportion to their masses. Strong viscosity gradients in proto- 
planetary disks offer a way to reconcile the circular orbits of Solar System gas 
giants with the eccentric orbits of currently known extrasolar planets. 

Subject headings: celestial mechanics — planetary systems — accretion, accretion 
disks 



1. INTRODUCTION 

Orbital eccentricities e, periods P, and semi-major axes a of extrasolar planetary systems 
are plotted in Figure 1. The rightmost four points at large P represent the gas and ice giants 
in our Solar System (Lodders & Fegler 1998). The leftmost cluster of points at small P reflect 
tidal interactions between planets and stars that erased whatever primordial eccentricities 
these systems possessed (Lin et al. 2000). In the intermediate range of periods, orbital 
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Fig. 1. — Orbital eccentricities and periods of 64 planets. The semi-major axis of the 
orbit is computed from the period using a central stellar mass of 1 Mq. TC represents the 
tidal circularization region, RRC the regime proposed to have witnessed repeated resonance 
crossings, and P the proposed outer passive region in which little or no migration occurred. 
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eccentricities can be strikingly large, typically exceeding those of giant planets in the Solar 
System by factors of ~2-20. 

Several theories have been proposed to explain these large eccentricities. Many en- 
counter difficulties when applied to the majority of systems, and one remains incompletely 
developed. A stellar binary companion can secularly drive a planet's eccentricity (Holman, 
Touma, & Tremaine 1997), but nearly all known extrasolar planets orbit solitary stars. Alter- 
natively, dynamical instabilities afflicting two planets formed at close separation on circular 
orbits can eject one planet while inducing a large orbital eccentricity in the remaining body 
(Rasio & Ford 1996; Weidenschilling & Marzari 1996). However, close encounters engineered 
in this fashion result also in planetary collisions, leaving a large proportion of planets on cir- 
cular orbits that are not observed (Ford, Havlickova, & Rasio 2001). In a third scenario, 
orbital migration of two planets during which the orbital period of one planet approaches 
that of the other ("convergent migration") can lead to resonant capture and eccentricity 
pumping. While convergent migration and resonance capture likely underpin the orbital ec- 
centricities of GJ 876b and GJ 876c, two planets observed to occupy a 2:1 resonance (Marcy 
et al. 2001b; Lee & Peale 2001), all other extrasolar planetary systems presently evince no 
mean-motion resonant behavior. A fourth explanation invokes gravitational interactions be- 
tween planets and the disks from which they formed (e.g., Artymowicz 1998). The present 
theory of sateUite-disk interactions has only been derived to lowest order in e (Goldreich & 
Tremaine 1980); it is inadequate when applied to extrasolar systems for which e's can be 
as large as 0.25-0.95. In the present theory, whether the disk damps or excites e depends 
sensitively on the distribution of disk material near the planet (Goldreich & Tremaine 1980; 
Papaloizou, Nelson, & Masset 2001). While this distribution is currently uncertain, Lee & 
Peale (2001) find that the disk must strongly damp eccentricities to reproduce the orbital 
parameters of GJ 876 and to avoid the fine-tuning problem of having the epoch of resonant 
capture coincide with the dissipation of the disk. 

This Letter proposes a fifth mechanism for exciting orbital eccentricities: repeated res- 
onance crossings of two planets migrating on divergent trajectories. In §2 we argue that 
divergent orbital migration of two planets is more likely than convergent migration. In §3 we 
demonstrate how divergent migration may lead to substantial eccentricity excitation. In §4 
we highlight the requirements and qualitative predictions of our theory and areas for future 
work. 
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2. Disk-Driven Divergent Drift 

Migration of planets can be driven by tidal interactions with their natal gaseous disks. 
Known extrasolar planets have sufficiently large masses (M > Af /, where Mj is the mass 
of Jupiter) that they clear annular gaps in disk material about their orbits (Ward 1997). A 
planet that opens a gap is thereafter slaved to the viscous evolution of its host disk, and 
undergoes so-called "Type II" drift (Ward 1997). The disk and its embedded planet at 
stellocentric distance r slide towards the star on the viscous diffusion timescale, 

2/ 2/ , c .r.Af r \i/2 /1000K\ /10-3\ 
tj, = t/\t\ ~ T^/y ~ T^/aCsh ~ 2 X 10^ j (^^^ j j ■ (1) 

Here v = aCgh is the viscosity of the disk, c^, /i, and T are the sound speed, vertical scale 
height, and temperature of disk gas, respectively, and < « < 1 measures the strength of 
angular momentum transport intrinsic to the disk. Equation (1) assumes a central stellar 
mass M* = Mq. 

The diffusion time Id cc ^jTa almost certainly increases with increasing r. Disk 
temperatures fall radially outwards as the stellar flux and the disk's absolute gravitational 
potential energy per unit mass diminish. Sources of angular momentum transport include 
(1) dissipation of density waves excited by numerous, densely nested planets that are insuf- 
ficiently massive (M < Mq, where is the mass of the Earth) to open gaps (Goodman 
& Rafikov 2001), and (2) the magnetorotational instabihty (MRI) that afflicts sufficiently 
ionized disks (e.g.. Stone et al. 2000). Mechanism 1 is capable of generating a. < 10~^, 
where the exact value depends on the spatial density of small, as yet undetectable planets. 
Mechanism 2 has been demonstrated to generate a ~ 10~^-10~^, the exact value correlating 
positively with the electrical conductivity of disk gas (Fleming, Stone, & Hawley 2000). The 
conductivity decreases with decreasing temperature, so that under mechanism 2, da/dr < 0. 

The standard MRI operates only in disk regions r < that are sufficiently hot, T > 
1000 K, that thermal ionization of trace metals and sublimation of dust grains permit the 
magnetic Reynolds number to exceed the threshold required for instability. Accretion disk 
models place between ~0.1 AU and ~1 AU (Gammie 1996; Bell et al. 1997; D'Alessio 
et al. 1998) — distances at which many extrasolar planets are presently located (see Figure 
1). In the absence of mechanism 1, it is possible that a is effectively zero at r > 1 AU. 
While Gammie (1996) has proposed that the standard MRI can still operate wherever gas 
densities are sufficiently low that Galactic cosmic rays can provide the requisite ionization 
levels, the likelihood of this prospect remains unclear for two reasons: (1) dust grains can 
severely reduce the electron density, and (2) even neglecting dust, and even if the magnetic 
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Reynolds number exceeds the critical value required for instability, the time required for a 
neutral molecule to collide with an ion is typically longer than a dynamical time, so that 
the bulk of the mostly neutral gas fails to accrete with the ions (Blaes & Balbus 1994). We 
return to the possibility of a static outer disk at the end of this Letter. 

What is important for what follows are the two recognitions that gap-opening planets 
at r < a few AU drift inwards at rates that arc (1) extremely slow compared to local orbital 
frequencies, and (2) different. Since dt£)/dr > 0, two gap-opening planets at r < a few AU 
drift inwards such that the ratio of the period of the outer planet, P2, to that of the inner 
planet. Pi, grows. The divergence of P2/P1 implies that a scries of mean-motion resonances 
will be crossed. During each resonance crossing, the orbital periods of the two bodies are 
momentarily commensurable; that is, the ratio of their orbital periods approaches and then 
exceeds a ratio of small, positive integers. 



3. Resonance Crossings on Divergent Orbits 

By contrast with the case where the period ratio P2/P1 converges towards unity, the 
probability of resonant capture is zero for diverging orbits (Sinclair 1972; Henrard & Lemaitre 
1983; Peale 1986; Yu & Tremaine 2001; and references therein). Nonetheless, as in the con- 
vergent case, each divergent passage can substantially alter the orbital eccentricities (Hen- 
rard & Lemaitre 1983; Peale 1986; Malhotra 1988; Dermott, Malhotra, & Murray 1988) and 
semi-major axes of the migrating bodies. 

We illustrate the underlying mechanics by considering the problem of a massive, in- 
wardly migrating planet about a star, and its effect on a massless test particle on a more 
distant, co-planar orbit. We take the mass of the planet Mi to equal 1.5 x 10~^M*, where 

= Mq is the mass of the central star. Referenced to a coordinate system fixed on the 
star, the initial orbital semi-major axis and eccentricity of the planet are 04 = 1 AU and 
Ci = 0, respectively. In addition to feeling the Newtonian force of gravity exerted by the 
star, the planet feels an additional drag force -F^rag = ^Miv/tdrag, where v is the instanta- 
neous velocity of the planet and tdrag = 1-6 x 10^ yr is the timescale over which ai decays to 
0. This prescribed drag force is introduced to simulate the effects of disk-induced migration 
and does not directly affect the planet's eccentricity (Papaloizou & Larwood 2000). Our 
value for tdrag is chosen to illuminate the evolution on timescales that are not too long com- 
pared to P2; larger values of tdrag are probably more realistic and will be considered later. 
The test particle is initially placed on an orbit having semi-major axis 02 = 1-35 AU and 
eccentricity 62 = 0, and is initially positioned at an angle A = 180° away from the angular 
position of the planet. The test particle feels only the gravitational attraction from the 
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Fig. 2. — Anatomy of resonance passages involving an inwardly migrating massive planet 

and a masslcss test particle on a more distant, co-planar orbit. In panel (a) the ratio of 
orbital periods of the test particle (P2) and of the planet (Pi) is plotted against time. Panel 
(b) plots the true anomaly of the test particle whenever the particle comes within 20° of the 
planet. Panels (c) and (d) plot the evolution of the test particle's eccentricity and semi-major 
axis. The resonances responsible for abrupt changes in the eccentricity are labelled. 
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planet and from the star. The subsequent positions and velocities of the planet and of the 
test particle are calculated using a variable-order, variable-step Adams numerical integration 
scheme (Hall & Watt 1976). Figure 2 displays the evolution. As the inner massive planet 
migrates towards the star, the mean eccentricity of the test particle undergoes 4 distinct 
changes. These changes occur when the period ratio P2/P1 equals 5/3, 2/1, 4/1, and 5/1. 
Concomitant changes in 02 occur at these epochs of resonance passage. At times t > 600 yr, 
interactions between the particle and the now-distant planet are negligibly small, and the 
particle is left on a more eccentric (e2 = 0.21) and shghtly expanded (02 = 1.47 AU) orbit 
compared to its initial one. 

These changes in 62 and 02 occur because during passage through a resonance, impulses 
of velocity are imparted to the test particle from the planet at specific phases in the test 
particle's orbit for extended periods of time. Accelerations felt by the test particle from the 
inner massive planet are strongest near times of conjunction, when positions of the star, 
planet, and test particle fall on a straight line and in that order. A p : q resonance for which 
the planet executes integer p circular orbits for every integer q elliptical orbits traced by 
the test particle is characterized by \p — q\ conjunctions which occur at \p — q\ values of the 
particle's true anomaly /2.^ True anomalies /2 near every conjunction are plotted in Figure 
2b. For example, during passage through the 2:1 resonance, conjunctions are repeatedly 
occurring at values of /2 concentrated in the interval between 270° and 360°. Conjunctions 
in this quadrant amplify 62 and 02 (Murray & Dermott 1999). Analytic expressions for 
the magnitudes of eccentricity jumps through first-order and second-order resonances are 
provided by Dermott, Malhotra, & Murray (1988). 

Accounting for the finite mass of the outer body does not change our conclusions qualita- 
tively. In Figure 3 we showcase 2 scenarios involving an inner body of mass Mi/M* — 2x 10~^ 
and an outer body of mass M2/M* = 1 x 10~^. To demonstrate that the substantial eccen- 
tricities that are excited in our simulation are caused by passages through resonances and 
not by the mere close proximity of these massive bodies, we do not impose any differential 
migration for the first 4 x 10^ yr of the integration. The osculating eccentricities of both 
bodies do not exceed 0.075 during this phase when 6r, the instantaneous distance between 
the two planets, can be as small as 0.400 AU. Only when the drag force is applied to the 
inner planet at t > 4 x lO^yr, using tdrag — 1 x 10* yr, do the orbits diverge. The eccen- 
tricities of outer and inner bodies then undergo resonant excitation to values of ~0.5 and 
~0.2, respectively, in the step-wise fashion characteristic of resonance crossings. The degree 
of amplification in 62 exhibited in Figure 3 is greater than that in Figure 2 primarily because 



^True anomaly measures the angular position of an object with respect to its periastron, and increases 
in the direction of the object's motion. 
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Fig. 3. — Resonance passages involving two massive planets. Masses of the inner and outer 
planets are Mi/M* = 2 x 10~^ and M2/M* = 1 x 10^^, respectively. The drag force is 
applied to the inner planet over tdrag = 1 x 10'' yr starting at tstart = 4 x 10^ yr. At t = 0, 
the ratio of semi-major axes is 02/01 = 1.5, the osculating eccentricities are both 0.005 (ab) 
or 0.03 (cd), and the planets are separated by an angle A = 180°. At t < t start: the planets 
mutually excite eccentricities of less than 0.075. Only after t > tgtart is differential migration 
introduced; the eccentricities of both objects become substantially excited through repeated 
resonance crossings. Panels (bd) also plot apastron distances Q — a{l + e) and periastron 
distances q — a{l — e). 
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the former is based on a longer tdrag- We have experimented with 20 different sets of initial 
orbital elements and find that eccentricities e > 0.2 are excited in all our test cases provided 
the initial P2/-P1 < 2 so that the 2:1 resonance is crossed. 

The eccentricity and semi-major axis jumps are due largely to resonant interaction 
and not to close encounters. The separation Sr attains minimum values of 0.325-0.400 AU 

(Figures Sab) and 0.306-0.400 AU (Figures 3cd) at t = 4000-5000 yr, whereas much of the 
amplification in 62 and 02 occurs over t ^ 5000-6000 yr when Sr > 0.400 AU. We have 
verified by a separate integration that freezing the divergent migration at t = 4356 yr in 
the case of Figures 3cd — thereby freezing 6r at its minimum value of 0.306 AU — does not 
generate further eccentricity or semi-major axis jumps at subsequent times, thus ruling out 
the significance of close encounters. 

4. Discussion 

We have established that divergent orbital migration of two planets can lead to signifi- 
cant eccentricity excitation without resonance capture. The effectiveness of this mechanism 
hinges on the detailed, and here umodelled, interaction between disk gas and the embed- 
ded planets. It is possible that the disk gas exerts torques on the planets on short enough 
timescales that the smooth navigation of the resonant separatrix is disrupted; conjunctions 
between the two planets may not occur at the same orbital phases for extended periods of 
time due to gravitational "noise" generated by disk gas. This is a concern that we defer to 
future hydrodynamic simulations of planet-disk interactions. 

While we have staged our scenario within a gaseous circumstellar disk, it seems possible 
that divergent migration may also be effected by planetesimal scatterings (e.g., Murray et 
al. 1998; Hahn & Malhotra 1999). In their simulations of the formation of the Oort Cloud and 
the Kuiper Belt, Hahn & Malhotra (1999) find that Uranus and Neptune divergently migrate 
and cross a 2:1 resonance (see their Figure 6), thereby undergoing eccentricity excitation. 
The excitation is limited, however, because the planetesimal masses which they employ are 
large enough that the divergence of the planetary orbits is not adiabatic. 

Our proposed mechanism operates most effectively when the initial orbits of the two 
bodies are sufficiently close that powerful resonances for which |p — = 1 (e.g., the 2:1 
resonance) are crossed during subsequent migration. Formation of planets at such proximity 
is not unreasonable; a single giant planet embedded within a circumstellar disk may induce 
the collapse of a second planet in the vicinity of the 2:1 resonance (Armitage & Hansen 1999; 
Bryden et al. 2000). 
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Disk-driven divergent migration of two gap-opening planets is most effective if a ring of 
viscous disk material remains present between the two bodies. Kley (2000) and Brydcn et 
al. (2000) have performed pioneering numerical simulations of two planets embedded in a 
disk resembling the minimum-mass solar nebula. While these calculations have tentatively 
shown that a ring can fail to be shepherded between two planets, so that planetary orbits 
converge rather than diverge, the outcome is model-dependent. The results of Bryden et 
al. (2000) suggest that if the mass of the ring is larger than the masses of the planets, or if 
the ring's intrinsic a is large so that planet-driven waves are efficiently dissipated near ring 
edges, then the ring can be confined; see, in particular, their model G. 

The resonance crossings mechanism for generating eccentricities demands that each 
planet have at least one other planetary companion, either in the past or today. If Mi < M2, 
then ei > 62, and vice versa. While the extrasolar planetary system v And violates this 
requirement (Butler et al. 1999; Chiang, Tabachnik, & Tremaine 2001), the system HD 
168433 (Marcy et al. 2001a) satisfies it. The absence of a second companion today in a given 
planetary system may indicate that the companion was accreted onto the star, either because 
its eccentricity grew too large by resonance crossings, or because "Type II" interactions drove 
the planet to its fiery destruction. A signature of planetary accretion would be an enhanced 
host star metallicity, but enhancements due to accretion of Jupiter-like giants during the 
first ~10^yr of the life of the star are too small to detect reliably (Murray et al. 2001). 

The condition that the planets initially share the same orbit plane may be relaxed. We 
would expect a non-zero initial mutual inclination, i, to be amplified in analogous manner to 
the way eccentricities are excited. Lifting the planet out of the plane of the disk may represent 
a means of survival against continued migration. However, the degree of amplification in 
i is expected to be less than that in e because inclination resonances are at least second- 
order {\p — q\ > 2) in strength (Murray & Dermott 1999). Calculating the relative orbital 
inclinations offers a means of testing our theory; future astrometric missions such as the Full- 
sky Astrometric Mapping Explorer (FAME) and the Space Interferometry Mission (SIM), 
in conjunction with stellar spectroscopic measurements, can place bounds on the degree of 
misalignment between orbital axes of eccentric planets and the spin axes of their parent 
stars. 

If the MRI is the sole source of viscosity late in the life of a protoplanetary disk, we would 
expect gap-opening, Jupiter-mass planets at distances outside a few AU to have suffered little 
to no migration within the primordial gas disk. While it is unconventional to think of T 
Tauri disks as having Q; = Oatr>rd~l AU, such a model does not appear to violate 
observation or theory. It would require the disk to contain ~0.01 Mq inside r < lAU, 
a condition for which the disk remains gravitationally stable. Timescales for accretion of 
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this much material could be as long as those observed, ~10^yr, if a ~ 10~^ at r ~ lAU. 
Given a static outer disk, the theory of eccentricity excitation proposed here would predict 
that orbits of giant planets at r > a few AU be nearly circular. Giant planet orbits in our 
Solar System conform to this expectation. We await the results of ongoing Doppler velocity 
searches and future space-based interferometric surveys for extrasolar planets to confirm 
whether the orbital architecture of the outer Solar System is indeed commonplace. 

We thank C. McKec for thoughtful comments on the manuscript and emboldening 
discussions, and P. Goldreich, R. Malhotra, G. Marcy, F. Marzari, N. Murray, E. Quataert, 
and S. Tremaine for encouraging and informative exchanges. 
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